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Abstract 

We discuss umbral calculus as a method of systematically discretiz- 
ing linear differential equations while preserving their point symmetries 
as well as generalized symmetries. The method is then applied to the 
Schrodinger equation in order to obtain a realization of nonrelativistic 
quantum mechanics in discrete space-time. In this approach a quan- 
tum system on a lattice has a symmetry algebra isomorphic to that 
of the continuous case. Moreover, systems that are integrable, super- 
integrable or exactly solvable preserve these properties in the discrete 
case. 

Key words: umbral calculus, difference equations, symmetries, integra- 
bility, quantum mechanics, discrete space-time. 
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1 Introduction 



A sizable literature exists on discrete quantum mechanics, that is on quan- 
tum mechanics in discrete space-time. We refer to a recent review for mo- 
tivation and for an extensive list of references [1]. There are many reasons 
for considering quantum systems in discrete space-time. One is that physi- 
cal space-time may indeed be discrete, involving an elementary length and 
some minimal time interval. Then continuous theories would only be ap- 
proximations to the real world. Another reason is the usual one: on a lattice 
one can avoid some of the divergence problems occuring in quantum field 
theories. On the other hand, some properties of quantum systems are lost 
in any discretization. The aim of this article is to discuss a discretization 
of space-time in which the Schrodinger equation is replaced by a difference 
equation. This is done in such a manner that many of the essential proper- 
ties of the continuous system are preserved. In particular, we preserve the 
group theoretical and integrability properties of the Schrodinger equation. 
This is true for the time-dependent, as well as the stationary equation. The 
free equations, as well as those with potentials, after discretization have 
symmetry groups, isomorphic to those of the continuous case. Lie point 
symmetries, after discretization, may however act at several points of the 
lattice. 

Another property that we wish to preserve is that of integrability, and 
also superintegrability. By integrability, for an n dimensional quantum sys- 
tem, we mean the existence of n well defined algebraically independent Her- 
mitian operators {X±, ...,X n } (including the Hamiltonian H) commuting 
pairwise. Superintegrability means that there exist further independent op- 
erators, {Yi,...,Yfc}, 1 < k < n — 1, commuting with the Hamiltonian H, 
but not necessarily with the other operators Aj, nor Yi [2]-[10]. 

Finally, we wish to preserve exact solvability in the discretization, i.e. 
the fact that for certain systems (like the harmonic oscillator, or hydrogen 
atom) it is possible to calculate all energy levels algebraically. 

A mathematical tool that we shall use for the study of symmetries and 
exact solutions of linear equations is the so-called " umbral calculus" . This 
calculus, which originated in 19th century with the work of Sylvester, Cayley 
and others, was used for a long time as a useful tool to derive combinato- 
rial identities (see, for instance, [11]). Nevertheless, it was only with Rota 
et al. [12]-[15] that this calculus was put on an axiomatic basis using the 
language of linear algebra of operators. In Ref. [16], the interested reader 
can find an up to date survey concerning the origins of umbral calculus 
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and its many applications in several branches of mathematics, like combi- 
natorics, functional analysis, algebraic topology, theory of special functions 
and orthogonal polynomials, etc. 

Umbral calculus has recently been used explicitly [17], or implicitly [18]- 
[21], to provide discrete representations of canonical commutation relations, 
specially in the context of exactly solvable and quasi exactly solvable quan- 
tum systems [22]-[24]. Linear differential equations have been discretized 
in a symmetry preserving manner using commuting difference operators 
[25]-[27]. An alternative approach [28] to symmetries of linear difference 
equations makes use of a discretized version of the prolongation theory of 
evolutionary vector fields. Finally, umbral calculus was used in an implicit 
manner, to obtain several different symmetry preserving discretizations of 
the linear heat equation [29]. 

Symmetries of difference equations, mainly nonlinear ones, have recently 
received a lot of attention (see e.g. [25]-[39] and references therein). What 
has emerged for purely difference equations is that in order to capture the 
essential features and usefulness of symmetries of differential equations it is 
necessary to make serious adjustments. Either one must go beyond point 
transformations to generalized ones [32]-[33], or one must use symmetry 
adapted and transforming lattices (as proposed initially by Dorodnitsyn) 
[34]-[39]. 

In this paper we follow the first approach. We consider a fixed lattice 
and use umbral calculus to obtain symmetries acting simultaneously on more 
than one point of the lattice. We apply this approach to quantum mechanics. 

2 Umbral Calculus 

To make this article self-contained let us sum up in Section 2.1 some known 
definitions and also some results proven as theorems by Rota [13], Roman 
[15] and Dimakis et al. [17]. 

We shall actually need umbral calculus on spaces of many variables in 
order to study multidimensional difference equations. However, for simplic- 
ity of exposition and notation we shall in this section restrict to the case of 
one variable x. 
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2.1 General Theory. 

Let ^ be the algebra of formal power series in a variable x, and V the algebra 
of polynomials in the same variable. The algebras will be considered over a 
field F of characteristic zero. This field in the subsequent considerations will 
be identified with KorC. An element of # is of the form 

= /(a). (2.1) 

k 

The operations defined in # are the addition of series 

oo oo oo 

^a k x k + ^b k x k = Y,( a k + h)x k , (2.2) 

k=0 k=0 k=0 

and the multiplication 

(oo \ / oo \ oo/m \ 

E * k x k E hx i = E E a ^-j * m - ( 2 - 3 ) 
k=0 / \l=0 / m=0 \j=0 J 

The algebra $ is also called the umbral algebra [15]. 

A polynomial sequence p n (x) G V is a sequence whose ra-th element is a 
polynomial of degree n. We will denote by C the algebra of linear operators 
acting on J or V. 

Definition 2.1 A shift operator T £ C is a linear operator such that 

Tp(x) =p(x + a), (2.4) 
where p(x) is a polynomial and a G F. 

Definition 2.2 An operator F G C is said to be shift-invariant if it com- 
mutes with all shift operators (i.e. with T for all values of a). 

Definition 2.3 An operator U is said to be a delta operator if it is shift- 
invariant and 

Ux = c^0, (2.5) 

where c G F. 

Using Definition 2.1 and the linearity of shift-invariant operators one 
can prove the following result. If U is a delta operator, for every c G F we 
have 

Uc = 0. (2.6) 



4 



Definition 2.4 A polynomial sequence p n (x), n = 0,1,2,..., is called a 
sequence of basic polynomials for the delta operator U if 

p (x) = l, Pn (0) = Vn>0, (2.7) 

and 

U p n (x) = np n ^i (x) . (2.8) 

It is easy to show that every delta operator has a unique sequence of basic 
polynomials. We will denote by 3 the one-to-one correspondence between 
basic sequences and delta operators. 

Let A be the algebra of shift-invariant operators, endowed with the usual 
operations of sum of two operators, product of a scalar with an operator, 
and product of two operators. We introduce a multiplication operation * : 
A x £ ^ £, defined by 

F * O = [F, O] = FO - OF, (2.9) 

where F is a shift-invariant operator and O G £. In particular, if x denotes 
the multiplication operator x : p(x) — > xp(x) , then F * x corresponds to 
what in the umbral literature is known as the Pincherle derivative of F. In 
this case we will write 

F' = F * x = [F,x] . (2.10) 
Using the * multiplication the Leibnitz rule becomes 

F*(fg) = (F*f)g + f(F*g) (2.11) 

and the Jacobi identity is expressed by 

F*G*H + G*H*F + H*F*G = 0. (2.12) 

Let us now consider a pair of shift-invariant operators: a delta operator 
U G £ and its conjugate operator G £, defined in such a way that the 
Heisenberg-Weyl algebra is satisfied: 

[U, xf3] = l. (2.13) 

It is possible to prove that, if U is a delta operator, then the inverse of 
U' exists ([13], p. 18). Therefore, the operator /3 is determined by the rela- 
tion: 

0={U')-\ (2.14) 
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To prove this it suffices to notice that 

1 = [U, x0\ = [U,x]0 = U'0, 

where the property [U,0\ = has been exploited. Eq.(2.14) follows. 

Let us present some specific examples of realizations of the conjugate 
operators U and (3 in terms of derivatives and shifts, respectively. 

Example 2.1. The continuous case. We have: 

U = d x , 0=1. (2.15) 

Example 2.2. The discrete case. The variable x is defined over an 
equally spaced lattice, with spacing a. Two of the most common choices for 
the discrete derivative are: 

a) The right discrete derivative. 

C/ = A+ = ^— ^, P = T~ 1 . (2.16) 
a 

b) The left discrete derivative. 

1 - T" 1 

U = A~ = , p = T. (2.17) 

a 

Other cases will be considered below. 

Using the * multiplication of eq. (2.9) it is easy to construct the basic 
sequence for the operator U. Let us introduce the polynomial sequence of 
operators 

P n =(x(3) n , n € N. (2.18) 
The delta operator U satisfies the relation 

[U, {x(3) n } = n (x/3)™- 1 , n G N. (2.19) 

This is an immediate consequence of the definition (2.13) and of the Leibnitz 
rule (2.11). A proof is obtained by induction. From (2.19) we immediately 
obtain: 

U*P n = nP n - 1 , n € N (2.20) 

This shows that {(x/3) n } nGN is the basic sequence for the operator U, under 
the * multiplication. 
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Definition 2.5 An umbral correspondence is a map 1Z : C — > C defined by 

(s/?i) n 4 (x(3 2 ) n , (2.21) 

where P\ = {{xj3\) n } and P% = {(x(3 2 ) n } are 6asic sequences of operators 
for two delta operators U\ and U 2 respectively. 

The umbral correspondence (2.21) naturally induces a correspondence 
between the two operators U\ and U 2 , according to the following scheme: 

(x^r ^ ( X f3 2 ) n 

3\ 3 1 (2.22) 

We shall also denote the induced correspondence between delta operators 
by the symbol 1Z. 

Systems of equations connected by the umbral map (2.21) share many 
algebraic properties. A particular case of the umbral correspondence is when 
U\ is the standard derivative d x , and U2 is a discrete derivative A. Then 
according to the scheme (2.22) we have 

(2.23) 



x n 




n 


n 


d x 


^ A 



(2.24) 

Let us observe that Definition 2.5 generalizes the notion of an umbral 
operator introduced in [13]: an umbral operator R : V — > V is an operator 
(in general not necessarily shift-invariant) which maps some basic sequence 
of polynomials p n (x) into another basic sequence q n (x) : 

Pn {x) & q n (x) • (2.25) 

Indeed, we observe that, since (3 is a function of shifts and any constant 
is invariant under the action of a shift operator, an umbral operator R is 
deduced from the action of 1Z simply applying the sequence of operators 
{xf3) n onto 1: 

(xfar ■ 1 4 ( X (3 2 ) n ■ 1. (2.26) 
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From (2.19) we also get 

Ui (x/3,) n ■ 1 = n (xftf -1 -1, i = 1, 2. (2.27) 

An important consequence is that the umbral correspondence (2.22) pre- 
serves commutation relations between operators in C. In particular, it pre- 
serves Lie algebras. 

Indeed, let A\ be a m-dimensional Lie algebra, generated by vector fields 

{ V i}i=l,...,m of the form 

v i = ^2a j (x 1 ,..,x p )d Xj . (2.28) 

3 

The umbral correspondence (2.22) maps A\ isomorphically into an algebra 
A2, generated by the vector fields {v^} i _ 1 m , with 

V Y = Y, a i (^^i-'^A) A *r ( 2 - 29 ) 
3 

This follows from the fact that the umbral correspondence (2.22) preserves 
the Heisenberg-Weyl algebra. 

2.2 Umbral Calculus and Linear Difference Operators 

The umbral approach reveals its power in the study of linear difference 
operators. 

For our purposes, namely the study of difference equations and their 
continuous limits, we shall need only two types of delta operators. The 
first is simply the derivative U = d x , with (5 = 1. The second is a general 
difference operator that has d x as its continuous limit. We put 

m 

U = A = - S"a k T*, l,meZ, I < m (2.30) 
k=l 

where au and a are constants and T a = Tis the shift operator of eq.(2.4). 
In order for A in (2.30) to be a delta operator, it must satisfy eq.(2.5). For 
any function / (x) G # eq. (2.30) implies 

k=l k=l 



8 



Using a Taylor expansion around a = we get 



1 00 f(?) M m 
A/ (x) = -V li^ V o fc A;9 (2.32) 

g=0 * k=l 
Choosing / (x) = xwe immediately see that eq. (2.5) implies 

m 

5> fc = 0, (2-33) 
k=l 

and YlT=i a kk = c. We require that in the continuous limit Abe the deriva- 
tive d x ; this implies c = 1, i.e. 



J2a k k = l. (2.34) 



k=i 

Eq. (2.30) involves m — I + 1 constants a^, subject to two conditions (2.33) 
and (2.34). To fix all constants a k we must impose m — l — I further condi- 
tions, for instance 

m 

lq = ^a k k q = 0, q = 2,3,...,m-l. (2.35) 
k=l 

Conditions (2.33) and (2.34) are necessary and sufficient for U = A to be a 
delta operator which has the derivative d x as its continuous limit. 

Definition 2.6 A difference operator of order p = m — l is a delta operator 
of the form (2.30) satisfying eqs. (2.33) and (2.34). 

Theorem 2.1 If the difference operator A of order (m — I) > 2 satisfies 

the supplementary conditions (2.35) it provides an approximation of order 
a m-l Q j ^ g d er i va ti ve Q x . 

Proof. We immediately have from eq. (2.32) and eqs. (2.33), (2.34) 
and (2.35) 

■m—l _ m _ 

A/ „~„ f W + l^TTTy. f{m -'- l> {x) S <^'"-'- K < 2 - 36 > 

v ; k=l 

QED 
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Remark. Formula (2.30) defines U as an operator parametrized by a, 
where a G F. It may happen that for specific values of a the operator U could 
involve less than p points, and consequently its order would be less than p. 
Once a representation of U as a difference operator is chosen, these points 
can be easily determined by solving a linear system of algebraic equations. 

By way of an example, let us consider the following equation [40]: 

A 3 / Or) + 3A 2 / Or) + Af (x) - f (x) = 0. (2.37) 

If A = 2^1, for a = 1 eq. (2.37) becomes 

/Or + 3)-/Or+l) = 0, 

which is of second order, in an appropriate domain. 

In the following, U will be assumed to be an operator of order p paramet- 
rically depending on a, and we shall omit the simple analysis of the specific 
cases in which the order could be less than maximal. 

Theorem 2.2 If A is a difference operator of order p, then A = T 3 A, j£Z 
is a difference operator of the same order. 

Proof. Let us first prove the result for j = 1. We have 

TA = -Y j a k T k+1 = - V a k T\ a k = a k _ v 



a L — ' a 

k=l k=l-l 



Hence 



m+l m 

k=i+i k=i 

m+l m m 

k ® k = ( k + ^ ak = kcik = 1 - 

k=l+l k=l k=l 

Thus, conditions (2.31) and (2.33) are satisfied for U and that is all that is 
needed. The proof for j = — 1 is analogous and for j arbitrary the result 
follows by induction. QED 

Conditions (2.35) are not shift invariant. However, once to and £are 
chosen equations (2.35) can always be imposed. Their solution depends on 
to and /, not only on the shift invariant difference m — l. 
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Theorem 2.3 The operator (3 conjugate to the difference operator A of eq. 
(2.30) is 



-i 



P=[J2 a kkT k ) . (2.38) 



,k=l 



Proof. Using eq. (2.14) we have 

f3=(A , y 1 = [A,x]- 1 . 

Moreover 

1 / m m \ m 

[A,x] = - \ y^a k {x + ka)T k -x^a fc T fe =Y j a k kT k 
a \k=l k=l J k=l 

and (2.38) follows. QED 

Examples of difference operators and the corresponding operators f3 are 
A + and A~ of eq. (2.16) and (2.17). Both are of order 1. Higher order 
examples are 

T^^ f T | T" 1 ^ \ 



A s = ^^, 0=[—L-—) , (2.39) 



A (ni) = _ J_ ( T 2 _ q T + 3 + 2T -i) ^ = (_ T " 3T r ^ 1 
Da v y \ 3 / 

(2.40) 

A (^) = --L^-ST+ST^-T- 2 ), (3 = ^_ T 2 -4T-4r- 1 + r- 2 y 1 

(2.41) 

The operators A s , A^ 111 ^ and A^ /y ) approximate the derivative to order a 2 , 
a 3 and a 4 respectively. 

Theorem 2.4 The expression 

P n (x) = (xp) n • 1 (2.42) 

is a well defined polynomial in x of order n with finite coefficients depending 
on a finite number of nonnegative powers of the shifts a for any difference 
operator A . The expression for P n is 

n 

P n (x) = A k a n ~ k x k , A n = 1, (2.43) 
k=i 
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where all coefficients A k are finite and depend only on the coefficients ak in 
the definition of A (see eq. (2.30)). In particular, they do not depend on a. 

Proof. Let us consider the difference operator A of eq. (2.30) and define 
the quantities 

m 

lj = J2a k k j , 70 = 0, 71 = 1, j = 0,1,2,.... (2.44) 

k=i 

Let us now prove eq. (2.43) by induction. Let P n be a basic sequence of 
polynomials for any A, as given by eq. (2.27). Thus we put 

n+l 

= ^TB a x a , (2.45) 

a=l 

and must prove that thejcoefficients B a are finite and depend on a in the 
proper way (i.e. B a = B a a n+1 ~ a , where B a is finite and does not depend 
on a). 

We rewrite eq. (2.27) (with n substituted by n + 1) as 

AP n+1 = (n + 1) P n - (2.46) 

The left hand side is 

m n+l n+l 

„k a— k— 1, 



^ m n+l n+l a , -. 

AP n+1 = -5>5> a (* + 6ar = 5> a £ I 

b=l a=l a=l k=0 ^ ' 



x cr 7 a _ fc 



n+l n+l , 
k=0 a=k ^ ' 

Comparing powers on the left and right hand side of eq. (2.46), we obtain 
a system of linear algebraic equations for the coefficient Bk : 

n+l , x 

E S « I ^ " fc_1 7a-* = (n + 1) (2.47) 



a=k 



The system (2.47) has a triangular structure. For = n + 1, we get the 
identity = 0. For k = n , only one term is present on the left and we get 
B n+ \ = 1 (since we have 70 = and A n = 1). The value k = n — 1 gives 

fl n + 1 



'7! 



a — K 1 - 2) • 
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In general, the system (2.47) implies 

n 

B n -j = <r j+1 » k A k (2.48) 

k=n— j— 1 

where the coefficients fik are easy to calculate, but are cumbersome (and of 
little interest), so we do not spell them out. QED 

Let us present the first few basic polynomials P k (x) = (x/3) fc 1 for arbi- 
trary A as given by eq. (2.30) with jj defined in terms of a k by eq. (2.44). 
We obtain 

P = (xpf -1 = 1 
p 1 = (x/3) 1 • 1 = x 
p 2 = (x/3) 2 . 1 = x 2 - a l2 x 
P 3 = (x/3) 3 • 1 = x 3 - 3a 72 x 2 - a 2 (73 - 3 7 f) x 
p 4 = ( x /3) 4 . 1 = x 4 - 6cT7 2 x 3 + a 2 (-4 73 + I57I) x 2 (2.49) 

+CT 3 (-74 + 107273 - 157 3 ) X. 

For a — > 0, we obviously reobtain the basic series (sequence) for A = d x . 

For A + = -^-we have only two values of aj, namely a\ = 1, ao = — 1, 
hence 7^ = 1, j = 2, 3, ... The polynomials (2.49) in this case reduce to the 
well-known factorial powers P n = x (x — a) (x — 2a) ... (x — (n — 1) a). 

2.3 Linear Difference Equations and Umbral Equations 

Let us introduce the notation f = f (x/3), i.e. to each function / (x) € $ we 
associate an operator / G £. We shall consider an operator equation of the 
form 

n 

Y^A k U k f = g (2.50) 

fc=0 

where U is a delta operator and (3 is its conjugate operator defined in eqs. 
(2.13)-(2.14). We assume that the operators A k and g can be expanded into 
formal power series in (x/3). 

Definition 2.7 An umbral equation of order n is an operator equation of 
the form (2.50) in which the operators A k and g are given. The unknown is 
the operator f . 
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If U is specified to be U = d x , then (3 = 1 and eq. (2.50) reduces to a 
differential equation of order n. If U is a difference operator, (2.50) is still 
an operator equation. Projecting both sides onto a space of functions, i.e. 
applying them to a constant, we obtain a difference equation. The order of 
the difference equations obtained projecting eq. (2.50) may vary depending 
on the structure of the operator A k (since it acts on the operator /) and 
on the choice of the operator A (and consequently of (3) in terms of shift 
operators. 

Let us first take U = d x , (3 = 1 in eq. (2.50). The obtained linear 
ODE will have n linearly independent solutions fi (x). We can expand them 
into formal power series about any point xo, where xq is not a singular 
point of the equation. Now let U = A be a difference operator and (3 the 
corresponding conjugate operator. Then /j {x(3) 1 viewed as a formal power 
series, will be a solution of the corresponding difference equation. 

Definition 2.8 We shall call f ■ 1 an umbral solution of the difference 
equation 

n 

Y / A k U k f-l = g-l (2.51) 

fc=0 

if the real valued function f (x) is a solution of the differential equation 

n 

Y,A k d k J{x)=g(x). (2.52) 

fc=0 

Thus each solution of the ODE (2.52) provides a formal power solution of 
the difference equation (2.51) (and of the umbral equation (2.50)). However, 
eq. (2.51) and (2.50) may have other solutions. Indeed, for a linear difference 
equation with constant coefficients we have the following theorem. 

Theorem 2.5 Let U be a difference operator of order p and let us assume 
that the operators A k in eq. (2.50) are constant. Equation (2.51) will then 
have np linearly independent solutions, n of them umbral ones. 

Proof. Eq. (2.51) in this case is a difference equation involving np + 
1 different points. Hence to obtain a solution in a new point we must specify 
initial conditions in points. This provides np linearly independent solu- 
tions, uniquely defined in the lattice points x n = xo+na [40, 41]. Now, let us 
consider the continuous limit of eq. (2.51). It is a linear partial differential 
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equation of order n, possessing analytic solutions which can be expanded 
around any nonsingular point. Applying the umbral correspondence to the 
series expansion of these solutions, we obtain n solutions of eq. (2.51) which 
are expressed as formal power series in (x(3) , and therefore are elements of 
the algebra 5- These are the umbral solutions admitted by eq. (2.51). The 
remaining (n — l)pdo not belong to QED 

When Ak are polynomials in (x{3), then additional shifts may appear in 
the explicit form of the equations coming from the umbral equation (2.50) 
via projection and their order may be different than np. 

As an example, let us consider the "umbral Airy equation" 

[A 2 + axj3] $ = 0, a = const. (2.53) 

For A+ = (T - 1) /a, (3 = T _1 and £ (x) = § ■ 1 we have 
1 



a 2 



(x + 2a) - 2 ^ (x + a) + (x) + ax ^ (x - a) = 0. 



This is a third order difference equation since it involves the function <]/ (x) 
at the points x + 2a, x + a, x and x — a. For A = A s eq. (2.53) would seem 
to involve infinitely many points: 

-l 



T + T~ 

^ (x + 2a) - 2 ^ (x) + ^ (x - 2a) + ax [ ) * (x) = 0. 



1 

(2.54) 

However, multiplying eq. (2.54) by [5 1 we obtain 



1 

4^2 



^ (x + 3a) — $ (x + a) — ^ (x — a) + ^ (x — 3a) + ax $ (x) = 

(2.55) 

This equation is a 6th order difference equation. 

As a simple example of umbral and nonumbral solutions, let us consider 
a first order homogeneous umbral equation with constant coefficients: 

Uf = af, a^O. (2.56) 

For U = d x , the solution is 

f(x) = Ae ax . (2.57) 

Now, let us consider the first order difference operator A + . Eq. (2.56) 
reduces to 

f(x + a)-f(x) = aaf(x). (2.58) 
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We look for a solution in the form / (x) = X x and find 

\ = {l + aa)°. (2.59) 
Thus we obtain a single solution 

= A(l + aa)° (2.60) 

and of course we have 

lim f(x) = Ae ax . (2.61) 

cr^O 

The umbral correspondence provides the solution 

f u (x) =Ae axT ~ 1 -1. (2.62) 

Expanding (2.60) and (2.62) in formal power series in a, we find that the 
two series coincide, i.e. f\ = f u . 

For comparison, let us consider the second order difference operator A s . 
Equation (2.56) in this case yields 

f(x + a)-f(x-a) = 2aaf (x) . (2.63) 

Putting / (x) = X x we obtain two values of A and the general solution of eq. 
(2.63) in this case is 

f = A 1 (Vl + a 2 a 2 + aa)~+^ 2 (-l)^ (\A + a 2 a 2 - aa)~ = A 1 f 1 +A 2 f 2 . 

(2.64) 

The first solution has e ax as its continuous limit. The second one does not 
have a limit for a — > 0. The umbral correspondence provides the solution 



f u (x)=Ae ax ( I± ^ 1 ) -1 



(2.65) 



(see eq. (2.39)). Expanding into formal power series in awe find f u = fi 
and /2 is nonumbral. 

The question arises whether an expression of the type 

f = e ax(S (2.66) 

is meaningful, at least in the sense of a formal power series. The problem 
is that for a general difference operator A, the expression for (3, given in 
eq. (2.38), is quite complicated. If we expand /3into a power series in T, it 
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will for m — I > 3 involve infinitely many shifts. Convergence problems may 
arise. Luckily, it is not eq. (2.66) itself that provides the umbral solution 
of a difference equation. Rather, it is the projection of the operator / onto 
a space of functions, or formal power series. The expressions that appear 
in the corresponding expansions are P n (x) = {xj3) n ■ 1 and these are finite 
polynomials in x, and in the shifts a, with well defined finite coefficients (see 
Theorem 2.4). As a matter of fact these are the basic polynomials for the 
difference operator A defined in eq. (2.30). 

It follows that if we know a solution of the umbral equation (2.50) for 
U = d x , and have 

~ f(n) ,q\ 
n=0 

then for U = A as in (2.30) the corresponding umbral solution will be 

* f W (0) 

n=0 

The matter of convergence is a separate issue. 

Umbral calculus and specially the umbral correspondence also provide 
us with a powerful tool with which to handle symmetries of linear difference 
equations, both ordinary and partial ones. On one hand, we can discretize 
a linear differential equation, in particular the linear Schrodinger equation, 
via the (multidimensional) umbral substitutions 

d Xi A Xi , Xi -5 Xif3 Xi , d t A t , t(3 t - 

Lie symmetries, both point and generalized ones, of linear differential equa- 
tions can be expressed in terms of commuting operators. Since the umbral 
correspondence preserves commutation relations, it will also preserve sym- 
metries. However, we may have more symmetries than in the continuous 
case due to the nonumbral solutions of the determining equations. 

3 Discretization of the Time Dependent Schrodinger 
Equation Preserving all Point Symmetries 

Before considering discrete space-time, let us first give a detailed and rig- 
orous analysis of the point symmetries of the time-dependent Schrodinger 
equation in continuous space-time. The results will be presented in a form 
well suited for the discretization. 
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3.1 Point Symmetries and Commuting Operators in Contin- 
uous Space-Time 

We write the Schrodinger equation in M n+1 as 

Lt/j = 0, 

1 n cP" 

L = id t -H, H = --J2g^ + V(^,t). (3.1) 

k=l k 

A local Lie point symmetry transformation is generated by a vector field 
that we write in evolutionary form [42] as 

v £ = Qdj, + Q*a r , (3.2) 

r> dtp dip 

Q = r] - T m-^ k - (3 - 3) 

The functions 77, rand depend on t, x , ^and ip* where the star denotes 
complex conjugation. These functions are determined from the requirement 

pr {2 \ E (LYO \w =L *^ = o= 0, pr {2 \ E (L>*) \ L ^ =L * r=0 = 0. (3.4) 

The following theorem will provide a basis for studying the symmetries 
of a nonrelativistic quantum system. 

Theorem 3.1 All Lie point symmetries of the time-dependent Schrodinger 
equation (3.1) are generated by evolutionary vector fields of the form (3.2) 
with 

Q = X (^,t)+iXiP, (3.5) 
X = i(r (t) 8 t + 6 t) <9* fc - ^ *)) , (3.6) 

£ k {~x, t) = 7^ x kT' - A k ixi + f k (t) , (3.7) 



1 



<P (^, t) = ^r"r 2 + x k f' k +g(t)+i [jt' - B\ , (3.8) 

where the prime denotes a derivative. The function \ {~z \t) satisfies the 
Schrodinger equation (3.1), A k \ = —A\ k and B are real constants. The real 
functions r(t), fk{t), g (t) and the constants A k i depend on the potential 
and satisfy the equation 

r (t) V t + & (1?, t) V Xk + t'V + X -T"'r 2 + x k f' k ' + g' = 0, (3.9) 
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with r 2 = J2k=i x 1- Moreover, the linear operator X commutes with L on 
the solutions of the Schrddinger equation 

[L,X]^ | LV , =0 =0. (3.10) 

Proof. Eq. (3.4) implies a system of determining equations. Those 
among them that come from terms involving derivatives of ip, e.g. ^x^xx ? 
tyxx, tpx, k > ldo not depend on the potential V(~x,t). From them we 
obtain the fact that the corresponding transformations are fiber preserving 
and linear (inhomogeneous). That is, Q has the form (3.3) with rand 
independent of V'and ip* ■ From the same equations we find that r depends 
only on t and that ( x ,t) are linear in x . Thus and </>have the form 
(3.7) and (3.8), respectively. 

Once these conditions are satisfied, only one determining equation re- 
mains, namely eq. (3.9), involving the potential in a crucial manner. 

To prove the commutativity relation (3.10) we use the compatibility of 
the two flows 

= % = Q, (Ml) 

where A is a group parameter and Q is the characteristic of the vector field 
(see eqs. (3.5)-(3.8)). Equating the cross derivatives iptx = ipxt and using 
the equation L\ {~x, t) = 0, we obtain 

[H,X\ij> = -iX t il>, (3.12) 

where i\) is any solution of the Schrodinger equation. This is equivalent to eq. 
(3.10). Simply stated: finding point symmetries of the Schrodinger equation 
is equivalent to finding linear selfadjoint operators X commuting with L on 
the solution set of L. QED 
Comments. 

1. For any potential (~x, t) the function x (~#\ i), the constant B and 
a constant g = go are solutions of eq. (3.9). Hence we always have a 
"trivial" symmetry algebra 

S(x)=x{~x',t)dTp + x*(~x,t)dTp*, L X = 

N = il>dj,+i/>*d r (3.13) 
E = i (V><9</, - i^>*d^*), 
due to the linearity of the Schrodinger equation. 
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2. Each symmetry generator v provides us with a flow that is by con- 
struction compatible with the time flow (3.1), that is, we can simul- 
taneously solve the equations (3.11). The fixed point dtp/dX = 
corresponds to group invariant solutions. 

3. While the result presented in Theorem (3.1) is quite simple and natu- 
ral, we have not found it explicitly in the literature, so we have sketched 
a proof. For other results on point symmetries of linear differential 
equations, see. e.g. [43], [44] and [45]. 

Let us consider the implications of Theorem 3.1 for special cases of the 
potential V(~x,t). We shall omit the operators (3.13) that are present for 
any potential V ( x , t). 

Let us first consider the free Schrodinger equation (for further discus- 
sions, see also [46]). For V(~x,t) = in 1R 3+1 the so called "Schrodinger 
group" was first obtained by Niederer [47]. For n arbitrary we obtain its 
generalization, i.e. the group Sch (n) . Its Lie algebra can be written as 

P = d t , D = 2td t + x k d Xk - - (ipdf + ip*dy* ) , 

C = t 2 d t + tx k d Xk - -t (Y><9^ + V*<V) + -J 1 " 2 ~ V>*<V) , 

L ik = Xid Xk - x k d Xi , P k = d Xk , (3.14) 
B k = td Xk + -x k (ipd^, - ^*<V) , 

E = i(i/}d^ - Y>*<V) . 
The Levi decomposition [48] of this algebra for n > 3 is 

£~[sI(2,R)®0(n)]9ff n (3.15) 

where the radical H n is the n-dimensional Heisenberg algebra. Explicitly 
we have 

sl(2,R)~{P Q ,D,C}, 0(n)~{L ik }, H n ~ {P k , B k ,E} . (3.16) 

We included the central element E explicitly in (3.14) since it appears in 
the derived algebra of the Schrodinger algebra (it is also present in (3.13)). 

Eq. (3.9) implies that for a general time independent potential V ( x ) 
we have only one additional symmetry generator to the set given by eqs. 
(3.13), namely time translations Pq = dt- 
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For a central potential V = V (r), the additional elements are time 
translations and rotations 

P = d t , L ik = xid Xk - x k d Xi , 1 < i < k < n. (3.17) 

In the case of a translationally invariant potential V = V (x n ), the ad- 
ditional symmetry elements are 

P = d t , Pj = d Xj , Bj = td Xj - ixj (ipd^ - ip*d^ ) , 1 < j < n - 1, 

L ik = Xid Xk - x k d Xi , 1 <i <k <n-l. (3.18) 

3.2 Symmetries of the discrete time dependent Schrodinger 
equation 

The umbral correspondence, together with Theorem 3.1 provide the tools 
necessary for a symmetry preserving discretization of quantum mechanics. 

Indeed, let us consider a discrete space-time, more precisely an n + 
1 dimensional orthogonal and equally spaced lattice with time step at and 
space steps a k , 1 < k < n. In this space we write a "Schrodinger difference 
equation" 

L d ^j = 0, L = iA t - H D (3.19) 
1 n 

Hd = --^2A XkXk +V (x 1 [3i,...,x n f3 n ,t(3 t ) 
k=i 

where we have 

[A Xk ,x k (3 k } = l, [A t ,tf3 t ] = l. (3.20) 

Each A Xk , A t is some chosen difference operator and (3 k , (3t are the corre- 
sponding conjugate operators satisfying the Heisenberg commutation rela- 
tions (3.20). 

The continuous limit of eq. (3.19) is eq. (3.1), obtained by taking 
a k ->■ 0, a t -> 0, i.e. 

A XkXk ^-^, At^Ot, ft-l, (3.21) 

Let us assume that in the continuous limit the obtained Schrodinger equation 
(3.1) is invariant under some Lie point symmetry group generated by some 
evolutionary vector field (3.2). Symmetries of linear difference equations on 
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fixed lattices can also be expressed in terms of evolutionary vector fields 
[28], [36]. For eq. (3.19) we put 

= Q D d^ + Q* D d r (3.22) 

Qd =VD- T D A t 1p - €k D A>x k ll> 

where rjo, td and ^k D are functions of Xifii, t(3 t , i> and ip* . The functions ip 
and 4>* are to be evaluated at the points Xifii, tfit- 

The prolongation of the vector field (3.22) must act on the dependent 
variables ip and -0* and on their discrete derivatives Afip, A XkXk ip. As in the 
continuous case, we require that an infinitesimal transformation 

XkPk = x k Pk, tfit = tfit 

4> {xk%M) = ^ (xkPk,tp t ) + XQ, A « 1 (3.23) 

should take a solution ip into a solution ipof the same equation (in new 
variables). First of all, we have 

fik = fik, fit = fit (3-24) 

since fik and fit are expressed in terms of shifts operators and we are con- 
sidering equations on a fixed (not transforming) lattice. Eq. (3.23) is an 
infinitesimal transformation in the evolutionary formalism, hence only the 
dependent variables transform. The transformation of the discrete deriva- 
tives is given by 

A t V = A t ip + AA t Q 

A XkXk $ = A XkXk i(j + XA XkXk Q (3.25) 

where A t , A Xk , etc. are discrete total derivatives. One can of course also 
introduce discrete partial derivatives [36], but we shall not need them here. 
In terms of the vector fields of eq. (3.22) the prolongation of is 

prvg = Q D d^ + Q l D d^ + Q x D kXk d AxkXk ip + ... + c.c. (3.26) 

where c.c. denotes the complex conjugate terms and we have 

Q l D = A t Q D , Q x ^ = A XkXk Q D . (3.27) 

The determining equations for the characteristic Q d are obtained as in the 
continuous case, i.e. from the invariance condition 

prvg (L D ip) \ LD i(,=L* D r=o= 0, P rv D ( L h^*) |l d v=l^*=° = °' ( 3 - 28 ) 
From this we conclude that the following theorem holds. 
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Theorem 3.2 The discrete time-dependent Schrddinger equation (3.19) al- 
lows a Lie algebra of "umbral symmetries" isomorphic to that of its contin- 
uous limit (3.1). This Lie algebra is realized by vector fields (3.22) with 



Qd = X (xkPk,t/3t) + iX D i/> 



X D = i 



T(tf3 t )A t + ^2t k A Xk 



1 n 
tk = -x k (3 k A t T - ^2 A kiXiPi + fk {tfh) 



(3.29) 
(3.30) 

(3.31) 



l=i 



-A tt T {xkPkf + x AA t f k + g (t(3 t ) 



k=i 



k=i 



+ i 



n 



- (A t r) - B 



(3.32) 

The function x sa tisfies the discrete Schrddinger equation (3.19), A k i = 
—Aik and B are real constants. The real functions r , f k and g all depend 
only on tfit and the potential V (x k (3 k ,t(3t) satisfies: 

n n n 

rAtV+Y, iikA Xk V+(A t T) V+- (A m r) £ (x fc /? fc ) 2 +^ x k p k A tt f k +A t g = 0. 



k=i 



k=i 



k=i 



(3.33) 

Finally, the difference operator Xp commutes with Lo on the solutions of 
the discrete Schrddinger equation (3.19): 



[L D ,X D ] ip \ L 



Dlp=0 



0. 



(3.34) 



Proof. The proof of Theorem 3.2 is quite analogous to that of Theorem 
3.1 in the continuous case. To see the similarities and differences, let us 
restrict ourselves to the case n = 2. 

The invariance condition (3.28) implies the following determining equa- 
tions 

€k,ip = £fc,v* = T i> = r V>* = 0> A Xi r = 0, </></,* = 0, = 0, (3.35) 

A X1 & + A X2 ^ = 0, (3.36) 
A t r - 2A xl £i = 0, A t r - 2A X2 & = 0, (3.37) 

XIX l 

6 + A 

X2X2 

Ci = o, 
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2iA t £ 2 + 2A X2 <f) 1 ^ + A X1X1 & + A X2X2 & = 0, (3.38) 

2ip {rA t V + ^A X1 V + &A X2 V + V A t r + V^} 

-2Vcp + 2iA^ + A XlXl (/> + A^^ = 0. (3.39) 

It is now obvious that eq. (3.30), (3.31) and (3.32) (for n = 2) provide a 
solution to eq. (3.35),..., (3.38) and that (3.39) reduces to eq. (3.33), once 
(3.35),..., (3.38) are solved. Eq. (3.34) then follows in exactly the same 
manner as in the continuous case. QED. 

Comment. 

There is an important difference between the continuous and the dis- 
crete case. In Theorem 3.1 we presented the most general solution of the 
determining equations. In Theorem 3.2 we presented a solution and added 
the requirement that the solution should have the correct continuous limit. 
Take for instance the function r. Equations (3.35) for any At, A x allow the 
solution r (t, fit), i.e. an arbitrary function of time t and the "shift" operator 
fit independently. For first order operators A ± (see Section 2), the function 
r will depend only on (t(3t). This follows from the determining equations, 
and agrees with the result of umbral correspondence. Moreover, the result 
will have the correct continuous limit. Thus the discretization and the con- 
tinuous equation have isomorphic symmetry algebras. For other choices of 
the discrete derivatives we may get more general solutions. For instance, let 
us consider the case of a "symmetric" derivative: 

^ T= r.-^' T= r(x + g )-r(x- a) 
x 2a 2a y ' 

Eq. (3.40) has the general solution 

T = T (t) + T 1 {t)e 12 tr, (3.41) 

x = x n = xq + na. (3.42) 

We see that eq. (3.40) actually allows an x-dependence in r. However, the 
second term in (3.42) does not have a continuous limit (for a = x n+ \ — x n —> 
0). 
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3.3 Examples 

As in the continuous case, for any discrete potential we have the "trivial" 
symmetries (3.13) (with x , t replaced by x k (3 k , tflt). 

Let us consider the case V = 0, i.e. a free quantum particle in discrete 
space time of dimension n + 1. Operators commuting with the operator L 
of eq. (3.19) are obtained from (3.14) by the umbral correspondence. We 
obtain the "discrete" Schrodinger algebra 

n 

P = A t , D = 2 (tfr) At + Y, (xkPk) A Xk ~ j W> + ^*d r ) 

k=i 

C = (tp t ) 2 A t + jr (tp t ) (x k f3 k ) A Xk - l - (tp t ) tyfy + rd r ) + 
k=l 

n 

-j Y ( x kPk) 2 (M/> - ip*dip*) 

k=l 

Uk = (xifii) A Xk - {x k (3 k ) A Xi , P k = A Xk , 
B k = (tfc) A Xk + % - (x k p k ) tyfy - Vd r ) 

E = i(^-rO r ). (3.43) 

For a general time independent potential V (xifti) we have only one ad- 
ditional (to (3.13)) symmetry generator, namely time translations Pq = A t . 

For a time independent central potential V = V (j2i ( x i(3i) 2 ^J the addi- 
tional symmetries are expressed by the operators 

P = A t , L ik = XiPiA Xk -x k (3 k A Xi , l<i<k<n. (3.44) 

For a translationally invariant potential V = V (x n /3 n ) the additional 
symmetry operators are 

= A t , 

Pj = A Xj , Bj = tp t A Xj - ixjpj -4>*A r ), 1 < j < n - 1 

Lik = Xi/3iA Xk - x k /3 k A Xi , l<i<k<n-l. (3.45) 
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4 Discrete Superintegrable Systems 



The umbral calculus provides a systematic method for transfering results 
from standard quantum mechanics to quantum mechanics in a discrete 
space-time. This is particularly simple if the results are formulated in 
terms of commuting differential operators. It has been shown elsewhere 
[5] that there is a direct relation between generalized symmetries in quan- 
tum mechanics and higher order differential operators, commuting with the 
Hamiltonian. Here we shall briefly sum up the results and then adapt them 
to the discrete case. 

4.1 Generalized Symmetries in Quantum Mechanics 

Let us consider the stationary Schrodinger equation in real two-dimensional 
Euclidean space 

H^ = E^, H = ~A + V(x,y), (4.1) 

and look for second order generalized symmetries in their evolutionary form 
V s (3.2) with characteristic Q satisfying 

Q = Q(X, V, ijj, l/j x , l()y, i) XX ,4> X y, tpyy) . (4.2) 

We require that the second prolongation of the vector field v E should anni- 
hilate eq. (4.1) on its solution space, i.e. 

pr (2) v E (H-E)ij)\ H ^ =E ^ = 0, pr {2 \ E (H-E)4>\ m=Eli> = 0. 

H4>*=E4>* H4>*=E4>* 

(4.3) 

If we also require that Q be energy independent, we obtain the following 
result. 

Theorem 4.1 The characteristic Q of the evolutionary vector field v E = 
Qd^p + Q*d^* , corresponding to a second order generalized symmetry of the 
Schrodinger equation (4-1) has the form 

Q = X^ + X (x,y) (4.4) 

X = aL 2 3 + b(L 3 P 1 + P 1 L 3 ) + c(L 3 P 2 + P 2 L 3 ) + d(P?-Pi) (4.5) 
+2ePiP 2 + aL 3 + j5P x + 7^2 + <P{ X , v)- 
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The function x (x , y) satisfies the Schrddinger equation (4-1) '• The operator 
X commutes with the Hamiltonian H 

[H,X}=0. (4.6) 

The quantities a, e, a, /?, 7 are constants and 

Pi=d x , P 2 = d y , L 3 = yd x -xd y (4.7) 

are generators of the Euclidean group E 2 . 

For a proof of Theorem 4.1, see Ref. [5]. 

The commutativity relation (4.6) is equivalent to the following linear par- 
tial differential equations satisfied by the potential V (x, y) and the function 

(f){x,y) 

[a (yd x - xd y ) + Pd x + 7 d y ] V (x, y) = 0, (4.8) 
(-axy -bx + cy + e)(V xx - V yy ) + [a(x 2 - y 2 ) - 2by - 2cx - 2d]V xy 

-3{ay + b)V x + 3(ax - c)V y = 0, (4.9) 
c/) x = -2(ay 2 + 2by + d)V x + 2{axy + bx-cy- e)V y , (4.10) 
4> y = 2{axy + bx-cy- e)V x + 2{-ax 2 + 2cx + d)V y . (4.11) 

Here eq. (4.9) is the compatibility condition for the two equations (4.10) 
and (4.11). Eq. (4.8) is easily solved. 

For a/0we can translate x and y to transform (3 — > 0, 7 — ► 0. Then 
the potential is rotationally invariant: V = V (r). 

For a = 0, /3 2 +7 2 / Owe can rotate to obtain (5 — > 0. Then the potential 
is translationally invariant: F = V (x). 

To avoid the geometric symmetries (4.7) we solve eq. (4.8) trivially by 
imposing a = (3 = 7 = 0. We then simplify the second order operator 
X of eq. (4.5) by rotations, translations and linear combinations with the 
Hamiltonian H. 

These transformations leave two expressions in the space of the coeffi- 
cients a, ...,e invariant, namely 

I 1 = a , I 2 =\{2ad-b 2 + c 2 ) 2 + 4(ae-bc) 2 ]. (4.12) 

In the nongeneric case when I\ = I 2 = 0, a third invariant exists, namely 

I 3 = d 2 + e 2 . (4.13) 
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Using these invariants, one obtains 4 equivalence classes of operators X and 
correspondingly, four classes of potentials allowing for the existence of an 
operator X, commuting with the Hamiltonian [3], [5]. The existence of 
one second order operator X , satisfying (4.6) makes the system integrable. 
Moreover, the corresponding Schrodinger equation will allow separation of 
variables in Cartesian, polar, parabolic or elliptic coordinates, which in the 
separable system depends on the values of the invariants (4.12) and (4.13). 

We are interested in the case of superintegrable Hamiltonians, when two 
operators X\ and X2 exist, satisfying 



[H,X 1 ] = [H,X 2 }=0, [X 1 ,X 2 ]^0. 



(4.14) 



Four classes of such potentials exist, each allowing the separation of variables 
in at least two coordinate systems. The Hamiltonians and corresponding 
integrals of motion are: 



1. 



Hi = ~\ ipl + ^) + ! V + V 2 ) + ^ + ^ (4-15) 



Xi = Pi - Pi - 



2/2 2\ 1 a b 



a b 

+ 



X 2 — /-,», 1 ., 1 . , 1 ■ 
\cos z <p sin z <f> J 

x = r cos (j), y = r sin < 



2. 



Hn = ~\ {dl + d^) + u\2J + ^) + ^ 2 +bx 



Xi = pi - pi 



u?{Ax 2 -y z ) +bx- 



y 



(4.16) 



X 2 = L 3 P 2 + P 2 L 3 - 2iv 2 xy 2 + - by 2 
The remaining two systems are best written in parabolic coordinates 

1 



x = - ( f - rj 2 ) , y = £77. 



(4.17) 



28 



3. 



Hiii 



1 1 



2 e 2 + r/ 2 v 4 ^ £ 2 + 7/ 2 V £ 2 



2a + A + ^) (4.18) 



X 2 = L 3 P 2 + P 2 ^3 + 



£ 2 + r/ 5 



(<? + 4)' 



(For 6 = c = 0, a^O this is the Coulomb atom). The system allows 
separation of variables in polar and parabolic coordinates (and also in 
elliptic coordinates). 



4. 



(a? + a?) + 



2a + b£ + en 



L 3 Pi + P ± L 3 + 



X 2 — -L3-P2 + ^"2-^3 + 2 



2( 2 + i) 2U T,/ £ 2 + t? 2 

6r? (r? 2 - e 2 ) + c£ (e 2 - V 2 ) ~ 4ar?£ 



(4.19) 



(e 2 + r/ 2 ) 

g(^ 2 -?? 2 )+? ? g (eg- 6??) 
(€ 2 + r?2) 



The equation separates in two mutually orthogonal parabolic coordi- 
nate systems, namely (4.17) and a similar system with x and y inter- 
changed. 

For o/O, 6 = c = we again obtain the Coulomb atom. 

We shall call the systems Hj and Hu the generalized isotropic and gen- 
eralized nonisotropic harmonic oscillators, respectively. Similarly, Hm and 
Hjv can both be called generalized Coulomb systems. 

The Schrodinger equations for Hm and Hjy can be rewritten as 

{ " \ (91 + d r)- E ie + ^) + ^ 2+ ^ = -a*, (4.20) 



_b_ 

2E 



(ri - 



2EJ 



1> 



(4.21) 
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9 b2 + c2 \l 

respectively. Thus, the system Hjjj is reduced to Hj with the energy (— E) 
and coupling constant to 2 interchanged. The system Hjy is reduced to a 
"shifted" harmonic oscillator. This interchange of the energy and a coupling 
constant has been called "metamorphosis of the coupling constant" [49]. 

4.2 Discrete Generalized Harmonic Oscillators 

The umbral correspondence immediately provides us with discrete versions 
of these systems. 

Let us first consider the potential V/of eq. (4.15). The discrete version 
of this system is: 



H? = ~\ (A* + Aj) + ^ \(x(3 x ) 2 + (yp y ) 2 } + £ (x/3,)" 2 + \ (y(3 y 



(4.22) 



with the integrals of motion 
X 1 = 



\^l+^ 2 (x/3 x ) 2 +a(x/3 x y 2 



--A 2 y + u; 2 (y(3 y ) 2 + b(yf3 y y 



and 



X 2 = {x(3 x A y - y(3 y A x ) 2 
a(l + (xp x y 2 (yP y f) + b(l + {x(i x f (y(3 y y 
Similarly, the discrete version of the system with potential Vu is 



(4.23) 



(4.24) 



H? I = - 1 -(A 2 x + A 2 y )+u J 2 



2(x(3 x ) 2 + -(yP y 



+ 2 W +bx(3 x . (4.25) 



The second order operators commuting with the Hamiltonian (4.25) are 



X l = Al - A 2 - 



(4 (x(3 x ) 2 - {y(i y f) + b {xf3 x ) - a (yP y y 2 ] , (4.26) 



and 



X 2 = [(yPy) A x - (x(i x ) A y ] A y +Ay [(yfiy) A x - (xfa) A y )-2u 2 {x(5 x ) (yf5 y f 



+2a{xf3 x ){yl3 y y 2 -b(y(3 y ) 2 . 



(4.27) 
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4.3 Discrete generalized Coulomb potentials 



To discretize the systems Hm and Hjy we again use the umbral correspon- 
dence, this time using parabolic coordinates. Thus, we replace 



(4.28) 



With these replacements it is a simple matter to write the discrete ver- 
sions of the systems corresponding to the potentials Vm&nd Vjv- Indeed, 
we have 



Hni 



l r 



A| + A 2 - 4a - 2b (£/%)" 2 - 2c (t/A,)" 2 ' 

(4.29) 



X 1 = [(^)A v -( v p v )A^} 2 -2 



(4.30) 

2 



X2 = [mf + W 2 ] 1 {( V f3,) 2 A| - (£/%) 2 A 2 + 2a [(£/%) 2 - (^ 

-26(r//3„) 2 (C/%)" 2 + 2c (^) 2 (r/^)~ 2 } (4.31) 



and 
Hjy = 



' + W 2 } 1 [ A I + A ' - 4a - 26 (£/%) - 2c (r/^)] 

(4.32) 



Xx = 



(£/%) 2 + W) 2 { (?/%) W (A| + A 2 ) + 



-6(r//3,)+c(^)] (ft^-foA,) 



4a 



(WW)} 



A| r (4.33) 



X 2 = 



(£/%) 2 + W) 2 ] 1 {(r/A,) 2 A 2 - (£/%) 2 A 2 + 2a [(£/%) 2 - (r^) 2 



(4.34) 



5 Exact Solvability and Spectral Properties of Dis- 
crete Superintegrable Systems 

We have shown that certain important properties of the Schrddinger equa- 
tion, such as point and generalized symmetries, and hence also integrability, 
are preserved when we pass from continuous to discrete space-time via an 
umbral correspondence. 
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Another important property of some quantum systems is their "exact 
solvability". This means that their Hamiltonian can be transformed into a 
block diagonal form with finite-dimensional blocks. In other words, their 
complete energy spectrum can be calculated algebraically. In more mathe- 
matical terms, we give the following definition. 

Definition 5.1 A quantum mechanical system with Hamiltonian H is called 
exactly solvable if its Hilbert space S of bound states consists of a flag of finite 
dimensional subspaces 

5 c5!cS 2 c ...S n c ... (5.1) 
preserved by the Hamiltonian 

HSi C Si. (5.2) 
All known exactly solvable systems also have the following properties. 

1. In appropriate coordinates and in an appropriate gauge, the bound 
state wave functions VPjv(~a?) are polynomials: 

*iv (~&) = 9 Pn {~s) , Si = SiCx). (5.3) 

The gauge factor g (1?) is a priori defined and can be energy depen- 
dent. The function Pjy ( s ) are polynomials of order N in the variables 
Sj. The integer N labels the subspaces Si in the flag. 

2. In the same gauge g and same variables Sj the Hamiltonian H can be 
written as 

H = ghg-\ hP N = E N P N (5.4) 

with 

h = a%T« k + ag m 2S2* (5.5) 

where a" fc and a^ lm are constants (subject to some further conditions 
[50]) and 

Tg = sfd Sk , a = 0,1, i,k = l,...,n. (5.6) 
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In other words, the gauge rotated Hamiltonian his an element of the 
enveloping algebra of an affme Lie algebra aff (n, R) (or one of its subalge- 
bras). It is clear that (5.5) guarantees that the Hamiltonian H will preserve, 
or decrease the order of the polynomials P/y. This is a concrete realization 
of the flag condition (5.2). 

We mention that all known quadratically superintegrable systems are 
exactly solvable, in particular those of Section 4 [6]. For the generalized 
harmonic oscillators the gauge factor g is equal to the ground state wave 
function and is energy independent. The generalized Coulomb systems have 
been reduced to the harmonic oscillators ones (see (4.20) and (4.21)). How- 
ever, due to the interchange of the energy and the coupling constant, the 
gauge factor g will be energy dependent. 

The aim of this Section is to show how exact solvability manifests itself in 
discrete space-time. First of all, let us consider an arbitrary one-dimensional 
linear spectral problem 

L(d x ,x)if>(x) = \il>(x). (5.7) 

Let x = be a regular point of this equation. Then any solution can be 
expanded into a Taylor series 

oo 

ij(x) = Y j a k x k . (5.8) 

k=0 

Using the umbral correspondence we write the umbral equation 

L (A, xf3) tp (x/3) = Xtp (x(3) (5.9) 

with the same eigenvalue A as in the ODE (5.7). Viewed as a difference 
equation, eq. (5.9) will have a formal power series solution 

oo 

V(x/?)1 = J>£(x/3) fc -l. (5-10) 

fc=0 

In particular, if (5.8) is a polynomial solution, then eq. (5.10) will also be a 
finite sum of terms involving the basic polynomials of the operator A. Thus, 
(5.10) will also be a polynomial and all convergence problems disappear. 

Now let us turn to the specific case of the generalized harmonic oscillator 
system with Hamiltonian Hj (see eq. 4.15) and its discretization (4.22). The 
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gauge factor g of eq. (5.3) and (5.4) is 
u (x 2 + y 2 ) 



g = x pi y P2 exp 



a = Pi (pi - 1) , b = p 2 (p 2 - 1) • 



We put ujx 2 = si, toy 2 = s 2 and in these variables the aff (2,R) operators 
of eq. (5.6) reduce to 

Ji = d Sl , J 2 = d S2 , J3 = sid Sl , J 4 = s 2 <9 S2 , 

Jb = S2d Sl , J6 = sid S2 . (5.12) 

The gauge rotated hamiltonian h and gauge rotated integrals of motion 
x~{ = gXig^ 1 , x~2 = gX2g~ 1 can now be written as [6] 

h = -2J Z J Y - 2J 4 J 2 + 2J 3 + 2J 4 - (2pi + 1) Ji - (2p 2 + 1) J 2 
= 2J 3 Ji - 2J 4 J 2 - 2J 3 + 2J 4 + (2pi + 1) Ji - (2p 2 + 1) J 2 
x5 = 4J 3 J 5 + 4J 4 J 6 - 8J3 J 4 + 2 (2pi + 1) J 5 - 2 (2p 2 + 1) J 3 + 

-2(2pi + l) J 4 + 2(2p 2 + l) J 6 . (5.13) 
By construction, all three of these operators will conserve the flag of poly- 
nomials 

Pn s 2 ) = ((sif 1 (s 2 f 2 \0<N 1 + N 2 <n) (5.14) 
and this is the reason why the superintegrable system with Hamiltonian Hi is 
exactly solvable. The actual solutions of eq. (5.4) are Laguerre polynomyals: 

HPnin = EnmPnrm Emn = Tl -\- 771, (5.15) 

Pnm (X,V) = 4-V2+P1) ( wx 2) (wy 2) . 

The umbral discretization will preserve the above properties and will 
give umbral Laguerre polynomial expressed in terms of x[3 x and y(3 y . The 
algebra aff (2, R) is represented by difference operators 

Ji = A sl , J 2 = A S2 , J 3 = (si/3i) A S1 , J 4 = (s 2 /3 2 ) A S2 , 

J 5 = (s 2 /3 2 ) A S1 , J 6 = (sift) A S2 . (5.16) 

The formulas (5.13) remain the same (with J« — > Jj) and all commutation 
relations are preserved, as are polynomial solutions. For similar results 
formulated in terms of operators acting in Fock spaces and the notion of 
isospectral discretization see Turbiner et al. ([19]-[24]), and [17] for further 
discussions. 
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6 Conclusions 



Much, if not all of nonrelativistic quantum mechanics can be viewed as the 
"theory of the enveloping algebra of the Heisenberg algebra". 

Indeed, let us define the Heisenberg algebra H n by the relations 

[X j ,Y k ]=6 jk C j,k = l,...,n (6.1) 

and then put 

Xj = Xj, Y k = -ihd Xk , C = ih. (6.2) 

We can say that all quantum mechanical operators lie in the enveloping 
algebra of H n , or in an extension of the enveloping algebra obtained by 
adding all formal power series in x±,...,x n , pi,...,p n . 

If we replace the coordinates Xj and the momenta pj by some other quan- 
tities satisfying the relations (6.1) then all polynomials and all power series 
in these objects will commute in the same way as the corresponding quantum 
mechanical quantities. 

Indeed, the umbral correspondence Xi — > Xifc, d Xi — > A Xi preserves the 
commutation relations (6.1) between quantum mechanical operators. Thus, 
the umbral correspondence allows us to consider quantum mechanics on a 
lattice and to preserve all properties of quantum mechanics in continuous 
space and time that are expressed in terms of the commutation properties 
of physical quantities (quantum mechanical operators). In particular, in- 
finitesimal point symmetries are preserved, as shown in Sections 3 and 4, 
respectively. Exact solvability is reduced to an algebraic property and then 
discretized in Section 5 (and also in the papers [20], [23]-[24]). 

Some nonalgebraic properties are lost in the discretization. For instance, 
the "umbral vector fields" (2.29), obtained by the umbral correspondence, 
do not generate global transformations (like rotations, or dilations). The 
solutions of umbral equations (obtained by the umbral correspondence) are 
often formal, i.e. they may diverge. 

The physical content of this article is based on the results contained 
in Section 2, where we presented and proved several theorems that, to our 
knowledge, extend the previously known umbral formalism to linear differ- 
ence equations. In particular we associate with a linear differential equation 
an abstract operator equation, written in terms of delta operators, using the 
umbral correspondence (2.23). Any representation of A and (3 in terms of 
shift-invariant operators provides a difference equation, whose analytic so- 
lutions can be obtained from the solutions / of the operator equation (2.50) 
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via the projection / • 1 = / (xj3) ■ 1. These are the umbral solutions admitted 
by a given difference equation. The other possible solutions do not have a 
continuous limit and are not provided by the umbral approach. 

Many of the results presented here can be considered also in the case of 
(/-difference operators [51]. The delta operator U in this case is defined in 
terms of a (/-shift operator satisfying 

T q f (x) = f (qx) . (6.3) 

In the simplest case, the (/-difference operator reads 

A, = j^-^ (T q - 1) , lim A q f = f x (6.4) 

and its conjugate operator, obtained imposing the Heisenberg commutation 
relation 

[A q ,x/3 q ] = l (6.5) 
is written in terms of a (/-shift and differential operators as 

(3 q = (q-l)(qT-l)- 1 xd x . (6.6) 

In this case we can still consider the umbral correspondence, but the 
A q operator is not a shift invariant operator, as [A q ,T q ] / 0. 

Among open questions, presently under consideration, we mention the 
following. 

• The umbral correspondence has lead us to various linear umbral equa- 
tions and difference equations. It would be of considerable interest to 
study their solutions directly and in the case of polynomial solutions, 
establish their relation to orthogonal polynomials of discrete variables, 
known in the literature [11], [52]. 

• Although many of the results presented in this paper can also be car- 
ried out in the (/-difference case [51], [53], [54], a reformulation of the 
umbral theory in this case is necessary. 

• The simultaneous diagonalization of commuting sets of second order 
differential operatores is intimately related to the separation of vari- 
ables in the Schrodinger equation. It would be important to investigate 
common solutions of commuting sets of difference operators from this 
point of view. 
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• A related problem is that of establishing a connection between umbral 
formalisms introduced in different coordinate systems, e.g. the carte- 
sian quantities A x , A y , (3 X , (3 y and the corresponding polar ones, or 
parabolic ones Ag, A^, /^introduced in Sections 4 and 5. 
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